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From a heuristic calculation of the leading order essential singularity in the distribution of Yang- 
Lee zeroes, we obtain new scaling relations near the ferromagnetic-Griffiths transition, including 
the prediction of a discontinuity on the analogue of the critical isotherm. We show that experi- 
mental data for the magnetization and heat capacity of Lao.rCao.sMnOs are consistent with these 
predictions, thus supporting its identification as a Griffiths ferromagnet. 
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The influence of disorder on ferromagnets remains, af- 
ter more than 30 years of effort, a complex and poorly- 
understood phenomenon. In its simplest form, disorder 
can be represented as a random spatial variation of the 
exchange interaction J in the bonds between neighbour- 
ing sites on a regular lattice. If a great enough fraction 
p > Pc of the bonds have J — 0, then one would expect 
that there is a vanishingly small probability of flnding 
a percolating pathway of bonds throughout the system, 
and the cooperative ferromagnetic phase would cease to 
exist. For smaller values of p, we would expect that the 
ferromagnetic phase will exist in a form weakened by the 
shortage of percolating paths; hence thermal fluctuations 
will destroy the ferromagnetic phase at a temperature Tc 
which is lower than the critical temperature Tq of the 
pure ferromagnet. However, as Grifflths showedQ it is 
not the case that the phase for Tc < T < Tq is purely 
paramagnetic, because in the thermodynamic limit, there 
can exist arbitrarily large volumes of the system that 
are devoid of disorder, with a probability exponentially 
sensitive to the volume. As a result, the free energy is 
non-analytic in external field, h, throughout the whole 
Griffiths phase. The effect of disorder is to partition the 
pure system into small ferromagnetic clusters. Depend- 
ing on the size, each cluster has a different value of Tc, 
so that the system as a whole exhibits a spectrum of Tc, 
spanning from the critical temperature of the pure sys- 
tem, Tg, due to arbitrarily large clusters, to some value 
of Tc, contributed by smaller clusters. 

Here we are concerned with the phase transition be- 
tween the ferromagnetic and Griffith's phases. Just as in 
the case of a pure ferromagnet, one would like to pre- 
dict the critical phenomena, but the non-analytic na- 
ture of the Griffith's phase makes it difficult to apply 
off-the-shelf renormalization group techniques 0, 0, Q 
or to posit simple scaling laws, despite recent theoret- 
ical progress S IE E|- Indeed, it is currently 
controversial whether or not there is clear experimental 
evidence O E E Q supporting the existence of the 
Griffiths phase. From the practical perspective, perhaps 
the most unsatisfactory aspect of efforts to relate the- 
ory to experiment is that critical exponents derived from 



conventional scaling laws are unrealistically large: for ex- 
ample, the critical isotherm exponent was rece ntlyli] 
estimated as (5 = 17. The breakdown of conventional 
scaling strongly suggests that the functional form of the 
scaling relations actually reflects the essential singulari- 
ties intrinsic to the Griffiths phase, and that some form 
of exponential scaling, rather than algebraic scaling, is 
appropriate. 

The purpose of this Letter is to address these problems 
by exploring the expected form of scaling relations that 
would follow from a leading essential singularity contri- 
bution to the statistical mechanics in the Griffiths phase. 
Such a contribution can be conveniently represented us- 
ing the Yang-Lee theory of phase transitions 0| to derive 
the scaling behavior in the Griffiths phase from a simple, 
physically-motivated ansatz for the distribution of par- 
tition function zeros, following arguments originally due 
to Bray and HuifangQ for the case of long-ranged ferro- 
magnets. We demonstrate that the leading singularities 
in the thermodynamics can be deduced, and that our pre- 
dictions consistently describe high quality magnetic and 
thermodynamic dataflll E| on the disordered Heisen- 
berg ferromagnet Lao.rCao.aMnOs. We emphasize that 
our purpose is only to identify the leading essential sin- 
gularities, and that it is beyond the scope of our work to 
provide a full description valid outside of the asymptotic 
critical regime. Nevertheless, the experimental data have 
sufficient resolution that we have strong support for our 
scaling predictions in this asymptotic regime. 

Magnetic properties are well-accounted for by our ap- 
proach, but heat capacity data are not expected to fol- 
low a simple scaling form-and indeed do not-due to the 
complex form of the theoretical predictions which arise 
from even the simplest Yang-Lee zero distribution func- 
tion that we use. The non-analyticity of magnetization 
in external field - a signature of the Griffiths phase - is 
also explicitly demonstrated. Our results provide strong 
evidence for a Griffiths singularity, and highlight the need 
for a more systematic renormalization group approach to 
understanding the singularities in such disordered sys- 
tems. 
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Yang-Lee zeroes and critical phenomena:- In 1952, Lee 
and Yang^^ developed a theory of phase transitions 
based upon the density g of zeroes of the grand partition 
function as a function of the complex fugacity and showed 
that the zeroes lie on a unit circle in the complex plane, 
parameterized below in terms of the angle 9. The distri- 
bution g(9) varies with temperature T and dictates the 
functional form of thermodynamics. Near a critical point, 
it is expected that g{9, T) exhibits behavior which reflects 
the non-analyticity of the thermodynamics^^) E|) and 
we provide this connection explicitly here, for both the 
case of conventional ferromagnetic critical point scaling, 
and then for the scaling near the Griffiths point. 

We begin with the scaling of the magnetization per 
spin, M(H, t), where H is the external magnetic field and 
t = [T — Tc)lTc is the reduced temperature, for a reg- 
ular Ising ferromagnet. The exact relationship between 
M{H,t) and g{9,t) can be written as 16] 



M{H,t) 



2n I d9 

10 



g{e,t) tanh(^i//fcsr)[l + cot2(6l/2)] 
1 + {ta.nh{fiH/kBT) cot(6'/2))2 



(1) 

where fi is the magnetic moment of individual spins. To 
extract the scaling behavior, we proceed by expanding 
around the critical point H = and t — 0. Because the 
singular behavior arises from the limit 9 —^ 0, we also 
expand the integrand about 9 = 0, resulting in 



M{H,t) = 2fj. 



d9 



(2) 



which is valid near the critical point, up to corrections 
reflecting a smooth background. A change of variables, 
9 = nHcp/ksT, gives 



M{H,t) =4fi 

Jo 



TT^^ 



where the upper limit is replaced by oo as _ff — > 0. This 
is the primary relation between M(H,t) and g{9,t) near 
the critical point. The scaling form of the magnetization. 



M{H,t) = \tffM{H/\tn, 



(4) 



where is an unknown scaling function, and (3 and 

6 are two critical exponents, implies a scaling form for 
g{9,t). Substituting Eqn. (gj) into Eqn. Q, Mellin trans- 
forming the expression and using the corresponding con- 
volution relation, we arrive at the scaling form of g{9,t) 
as 6* ^ 



ji9,t) = \tfG{9/\tn 



(5) 



where G{x) is a scaling function for g{9,t). 

Apart from exhibiting the scaling form of g{9,t) for 
normal ferromagnets, this exercise also shows that knowl- 
edge of g{9, t) can be gained by studying the scaling form 
of M{H,t), and vice versa. In the following section, we 
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FIG. 1: Data collapse of magnetization, M{h,t). Dots are 
experimental data and the line is the theoretical prediction 
for the universal scaling function, with go = 0.5, Aq = 0.5 
and /3r = 0.8. Data shown are in the range of t < 0.35 
and h = jiH/kBTc < 0.39. The insert shows the data on 
logarithmic scales. 



study the scaling behavior of M(H,t) in the Griffiths 
phase using a heuristically-derived form for g{9). We will 
see that the result is different from that of the pure case, 
reflecting the intrinsic essential singularity that charac- 
terizes Griffiths phases. 

Density of zeros for a disordered ferromagnet:- We start 
with the scaling form of g{9,t) derived on the basis of 
heuristic arguments by Bray and Huifang^ for disor- 
dered ferromagnets with short-ranged interactions: 

^ oo 

g(9, t) = -ny^ cxp(-A(t)r) tanh \r(i9 + e)] (6) 

TT ^ 

r— 1 

where e ^ 0+ and A{T) - (T - T^)^-'^'- as T ^ and 
A{T) ^ oo as T — > Ta- The exponent (3r is the order 
parameter exponent for the random case, and its value 
will reflect the universality class of the magnet, be it 
Ising, Hciscnberg, 0(n) etc. In the limit that 9^0 and 
A — > 0, the summand is dominated by peaks at values of 
r = (2n -f- 1)tt/29 for all non-negative values of n, whose 
height is given by {29/e{2n + 1)tt) exp(-(2n -\- 1)ttA/29). 
The width of these peaks therefore scales in the same 
way as their separation, both being proportional to 1/9. 
Thus, the expression takes the form: 



g{9,t)=goexp{-A{t)/\9\), 



(7) 



where go is a constant. The essential singularity in 
Eqn. ^ reflects the Griffiths phase of disordered fer- 
romagnets, not present in the pure case. Disordered 
magnets with long-range interactions have a power-law 
prefactor to the essential singularity but this is not 
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present in the short-range case. A disordered ferromag- 
net can be thought of as an ensemble of weakly interact- 
ing, finite-sized ferromagnetic clusters. When T ~ , 
only large clusters contribute to the overall magnetiza- 
tion. For each large cluster of linear size L, the small- 
est Yang-Lee zero is of the order of 6* ^ 1/mL'^, where 
m ~ [Tq—TY is the magnetization per spin of the cluster 
and d is the spatial dimension of the system. The prob- 
ability of a spin belonging to a cluster of size L follows 
the Poisson distribution, i.e., prob ~ exp(— cL''). As a 
result, large clusters contribute to g{0,t) in the form of 
Eqn. 0, where A{t) - (Tg - T)-/^ and /3 is the usual 
exponent for pure ferromagnets. 

For T < Tc, the system is in its ferromagnetic phase 
with nonzero magnetization, implying a nonzero value of 
g{e = 0,t)^. This requires A(t) = at T = to 
counteract the essential singularity at 6* = 0. For T ^ 
, we can expand A[t) and approximate it by A{t) ^ 
l^-Pr _ This accounts for the asymptotic behavior of A[t). 

With this form of g{6,t), the scaling behavior of 
M{h, t) can be obtained by substituting Eqn. Q into 
Eqn. ^ and making a change of variables, y — A{t)/d, 
yielding 



Mih.t) ^ gpAjt) 
a 2h 



, exp(-y) , . 

^,y/yy^ + {A{t)/2hY^ 



where we defined h = ^H/ksTc. Eqn. © can be 
written in terms of the exponential integral, Ei{x) = 
/tdt, as 



M{h, t) 



= -go'^[eyiY>{iA{t)/2h)Ei{A{t)/'K + iA{t)/2h)]. 

(9) 

We expand Eqn. Q about yl = 0, because A{t) = 
Aot^^^'' is asymptotically small in the critical region, re- 
sulting in 

^^ihA = -goe-^(*)/^x3[e^^(*)/2h^^(^^(^)/2/i)]. (10) 

This implies an approximate scaling form 
M{h,t) 



= exp(-A(t)/7r) [fM{A{t)/h) + 0{h)] , (11) 



where 



/m(x) = -gQ'^[e^^{ix/2)Ei{ix/2)], 



(12) 



and the corrections of 0{h) involve exponential integrals 
of A/h. This scaling prediction is valid in the limits A — > 
0, ft- ^ 0, but the ratio A/h has not been fixed by the 
analysis so far. 

Analysis of magnetization data:- We now analyze the ex- 
perimental data on Lao.rCao.sMnOa^, 0|, using the 
above results, to see if the data are consistent with the 
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FIG. 2: Experimental verification 
of the magnetization for A/h 
ear dependence of M/h on exp(- 
predicted by Eqn. 1131 . The ranj 
0.05 < A < 0.155. Legend: • h 
h = 0.0017, A ft = 0.0023, v h 
h = 0.0045, a h = 0.0057, o h -- 
h = 0.0122, -k h = 0.0141 and the 
prediction of Eqn. 1131 . 



of the asymptotic behavior 
oo. The data show a lin- 
A/tt)/A with slope 0.93 as 
16 of values of A plotted is 
= 0.0006, X h ^ 0.0011, o 
= 0.0028, + h = 0.0034, * 
= 0.0072, < h = 0.0090, > 
solid line shows the linear 



presence of a Griffiths phase. Fig. ^ shows the pre- 
dicted data collapse of the magnetization with fitted val- 
ues Tc = 218K, go = 0.5 ± 0.05, = 0.5 ± 0.05 and 
(3r = 0.8 ± 0.05. The error bars were obtained by esti- 
mating the best fit visually. The figure also shows the 
agreement between the theoretically-predicted universal 
scaling function and the collapsed data. Overall, the data 
scale quite well, and the scaling function of the collapsed 
data are close to that of the theory, except near the turn- 
ing point of the curve, where corrections to the leading 
order ansatz we have used for Eqn. (|ll|l become impor- 
tant, and the data are no longer in the asymptotic limits 
A and /i — > 0. In order to show this clearly, we 
calculate the asymptotics of the scaling function, in the 
limits A/h and A/h — > oo, where the leading terms 
in the magnetization can be calculated systematically. 



Asymptotic behavior:- From Eqn. (j^, 
havior of M {h, t), in the limit A/h ^ 



the asymptotic be- 
oo, is given by 



M{h,t) 2h 



-A/tt 



M50 



A 



1 - 



4ft2 



8h^ 
ttA 



O 



J 



(13) 

showing that the small field susceptibility, M/h, in the 
limit A/h — > oo, depends linearly on exp{— A/ tt) / A. This 
prediction recovers Curie- Weiss-like behavior, as verified 
by the experimental data shown in Fig. |(5J). The exper- 
imental data shows a slope of 0.93, consistent with the 
fitted range of values for 0.9 < 2go < 1.1. 

A profound difference between conventional and Grif- 
fiths ferromagnets is found in the limit A/h 0. It 



4 



is well understood, for conventional ferromagnets, that 
M{h,t) ~ h^^^ —I- as h ^ 0; this is, however, not true 
in Griffiths ferromagnets, as we can see by calculating 
the asymptotic behavior of M{h, t): 



M {h, t) IX 



In 



+ 7-1 



A 

2h' 



2h 



—+ —+0 (— 

(14) 

as A/h 0, where 7 0.57722 ... is the Euler- 
Mascheroni constant. This implies that M{h,t) 
figoTx/2 as h, A/h — > 0, i.e., a discontinuity of M{h,t) 
at h — on the critical isotherm. Experimental support 
for this surprising result is present in Fig. 3 of Ref. ([15j)j 
which documents the increase of the exponent S fitted to 
data assuming the conventional scaling applies. As the 
transition temperature is lowered by increasing disorder, 
the inferred value of S rises to as much as 16.9, incon- 
sistent with any known universality class, but indicative 
of a very rapid and dramatic rise in magnetization. In 
order to test the precise predictions made here, we show 
in Fig. © the experimental verification of Eqn. (|14|l . 
where the function M/figo + 2/i/7r — Ah/n^ is plotted 
against A/2h. The data points satisfying the criteria 
h < 0.0072 and A/h < 0.3 are shown in Fig. All 
the parameters were determined previously according to 
the data collapse in Fig. so that we have not made 
any additional fitting. The experimental data approach 
the theoretical curve as — > 0, and moreover tend to 
the universal number 7r/2 as dictated by Eqn. ifTH) . We 
conclude that the data are consistent with the predic- 
tion that M{h,t) is discontinuous at /i = in the limit 
A/h ^ 1, a prediction which follows from the essential 
singularity characterizing the Griffiths transitional. 

Heat Capacity:- We conclude with a brief discussion of 
the heat capacity, C{h,t). We can integrate Eqn. (|11|1 to 
obtain the free energy, F{h, t) and thence an expression 
for C{h, t). There is no prediction of data collapse, due to 
the interference from the exponential terms, in agreement 
with our failure to obtain data collapse from the data. 
This form contains singular terms of the form 



de e'^^^ ln(4/i2 9^) - exp(-A//i) (15) 



where the integral is restricted to the neighbourhood of 
the origin where Eqn. Q is valid and the upper limit 
is assumed to scale with /i, leading to the estimate of 
the essential singularity. Similar terms have also been 
predicted in Ref. (113), but more than the leading term 
must be retained in order to consistently compute the 
magnetization. 

In conclusion, we have argued that the essential sin- 
gularity of the Griffiths phase leads to novel features in 
the critical behavior, including a discontinuity of magne- 
tization in external field. These features are reproduced 
to the accuracy expected from our lowest order theoret- 
ical predictions by high quality experimental data from 
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FIG. 3: Asymptotic behavior of the magnetization for A/h 
0, showing that the data follow Eq. I14II . and exhibit a dis- 
continuity M figoiT/2 in the limit h 0. The data shown 
here are within the limit A/h < 0.3 and h < 0.0072. Legend: 
\/ h = 0.0028, o h = 0.0034, * h ^ 0.0045, + h ^ 0.0057, ■ 
h = 0.0072, X ft = 0.0090, a h = 0.0112, o ft = 0.0141. Solid 
line: theoretical prediction of Eq. I14II . 



Lao.7Cao.3Mn03, and lead to a consistent description of 
its critical behavior, supporting the identification of this 
material as a Griffiths fcrromagnet. 
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